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Thermal components in the early X-ray afterglow of GRBs 
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Abstract. The possible presence of thermal components in the early X-ray afterglows of 7-ray bursts is inves- 
tigated. We discuss both the presence of a thermal continuum and, in particular, of coUisional X-ray emission 
lines. We compute the predicted luminosity by a thin plasma for a range of metallicities for the continuum and 
the Ka lines of the elements Mg, Si, S, Ar, Ca and Fe. We show that light travel effects are dominant in the 
determination of the thermal continuum and line luminosities, and derive the relevant equations. We conclude 
that thermal lines and continua are unlikely to dominate the early afterglow of GRBs, unless the explosion site is 
surrounded by a very massive and extremely clumped shell of material. Such conditions are difficult to envisage 
in the close environment of GRB progenitor, unless they are excited by some strong precursor activity, like in the 
Supranova scenario. 
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1. Introduction 

Almost all the 7-ray bursts (GRBs) detected so far are 
associated to a transient X-ray afterglow (Lazzati et al. 
2002a). This afterglow is supposed to be due to the early 
deceleration of the fireball by the interstellar medium 
(Meszaros & Rees 1997) and its radiation produced by 
non-thermal synchrotron. Such an interpretation is cor- 
roborated by spectral (van Paradijs et al. 2002) and po- 
larimetric (Covino et al. 1999) observations. 

More recently, several X-ray afterglows observed by 
Newton-XMM had a spectrum that can be better fit by 
an optically thin thermal bremsstrahlung model (Watson 
et al. 2002) rather than with an absorbed power-law (even 
though the latter model cannot be unambiguously ruled 
out on purely statistical grounds). In particular, several 
high ionization emission lines were detected in the early af- 
terglow of GRB 011211 (Reeves et al. 2002, hereafter R02) 
which can be accounted for, together with the observed 
continuum, by a moderately enriched thermal plasma. 

The possibility of a thermal origin of the X-ray lines de- 
tected in several afterglows was first discussed in Lazzati et 
al. (1999; see also Vietri et al. 1999 and Kumar & Narayan 
2002), who pointed out how a 10 times solar plasma may 
produce a Fe Ka line with the observed large luminosity 
and equivalent width. In this paper we derive more rigor- 
ously the conditions for the environment and the heating 
mechanism that must be satisfied in order to observe line 
and continuum emission with the prescribed luminosities, 
equivalent widths (EWs) and variability time scales. 



2. Continuum and line emission in thermal plasma 

To compute the continuum emission from a thermal opti- 
cally thin plasma we adopt the treatment and basic equa- 
tions of Rybicki & Lightman (1979). To compute the lu- 
minosities of Ka lines from a thermal plasma, instead, we 
adopt the widely used code MEKAL (Mewe et al. 1985; 
Liedahl et al. 1995) as implemented in XSPEC (Arnaud 
1996). We concentrate in particular the elements Mg, Si, 
S, Ar, Ca and Fe, for which emission lines have been de- 
tected in the afterglow of GRBs. In Fig. |l| we show the 
line production efficiency for an optically thin thermal 
plasma as a function of the plasma temperature and for 
several values of metallicity. We plot Tyunc, i-e. the ratio of 
the Ka emission line for the six elements above, irrespec- 
tive of their ionization state, over the total luminosity of 
the plasma. Gray shading highlights regions in which the 
equivalent width (EW) of the lines is less than 100 eV, a 
robust lower limit to any emission feature detected in the 
afterglow so far. 

When using the MEKAL code to evaluate line emission 
luminosities, one has to remember that the code does not 
include any radiation transfer, since it assumes that the 
medium is optically thin to radiation. The actual optical 
depth of a cloud of plasma at temperature T depends on 
the temperature and on the frequency of the radiation con- 
sidered. When X-ray continuum radiation is considered, 
the optically thin approximation can be used up to col- 
umn densities Nh ^ <Tj^^ ~ 1.5 x 10^^ cm~^. If, however, 
line emission is concerned, it must be taken into account 
that intermediate-high Z elements retain some electrons 
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Fig. 1. Efficiency of line emission for several elements as a function of the temperature in an optically thin plasma. 
The parameter ?7iinc is defined as the line luminosity over the bolometric luminosity of the thermal plasma. Each panel 
shows 5 different lines. From top to bottom these refer to a metallicity 10, 3, 1, 0.3 and 0.1 solar. The gray shading 
highlights regions in the parameter space where an equivalent width lower than 100 eV is predicted. 



which may cause the plasma to be optically thick due to 
photoionization. In Fig. ^ we show the optical depth of a 
solar metallicity plasma with Nh = 1-5 x 10^^ cm~^ as 
a function of frequency for a range temperatures between 
10^ and 10^ K. The opacity of a cold gas is also shown for 
comparison. In the region of the considered emission lines, 
the plasma can be optically thick up to temperatures of 
several keV. This will limit the maximum line luminosity 
and EW: increasing the column density of the plasma will 
have no effect on the line luminosity since line photons 
will be able to escape freely only out of the optically thin 
surface layer of the medium. Instead of properly introduc- 
ing radiation transfer in the optically thin MEKAL code, 
which is a formidable task, in the following we will assume 
that, after the plasma becomes optically thick to radiation 
at the line frequency, the line luminosity does not increase 
if the column density of the gas is increased (see Lazzati 
et al. 2002b for discussions on a similar assumptions in 
reflection models). 

3. Thermal models for GRB afterglows 

Consider a cloud of plasma of electron density n at a tem- 
perature T which covers a solid angle — 2ti{\ — cos 6') 
along the line of sight to the observer at a distance R from 



the burst explosion site (see Fig. |^). The luminosity that is 
inferred from infinity, taking into account the light travel 
time effects, is given by|^ 

i « 1.7 X 10-23 T,i/2 EI — ^'^""^ + (1) 

* i?/c(l-COs6l)+thcat+tcool 

where EI is the emission integral {EI = neTii dV) , 
^cooi ^ 1-2 X 10^^ Tg^^2 n"-'^ s is the plasma cooling time 
and ^hcat is the heating time, i.e. the time during which 
heat (or energy) is supplied to the emitting plasma. It is 
assumed that during this heating time the temperature 
is held constant, i.e. the plasma is in equilibrium. In the 
following we will consider a uniform density plasma and 
approximate k, m ~ n. In this case EI — -n? V, where 
V is the volume of the emitting cloud. 

In many astrophysically relevant situations (think e.g. 
to the IGM in galaxy clusters) R/c ^ tcooi and both the 
heating and cooling time scales are irrelevant and Eq. |^ 
simplifies to the usual free-free equation: 

i = 1.7 X lO-^^Tg^^^;/ (2) 

^ Here and in the following we indicate a quantity Q as Q — 
W Q^. We also adopt CGS units. 
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Fig. 2. Photoionization optical depth of a = 1 plasma 
in thermal equihbrium as a function of the temperature 
and frequency. The range of frequencies of the considered 
Ka lines is shown by gray shading. 

In GRBs this may not be the case. In order to repro- 
duce afterglow observations in the non-cooling regime, 
the plasma cloud must satisfy simultaneously three con- 
ditions. First, it must produce the observed luminosity in 
the non-cooling regime 



erg 



(3) 



L = 1.7 X 10-23 Ts^/^n2^i?3 il 
^ 3 47r 

where the numeric value is typical for the early X-ray af- 
terglow (see e.g. GRB 011211, R02) and rm is the vol- 
ume filling factor of a possible shell-like or clumpy cloud. 
Second, it must fulfill the non-cooling condition: 



i?(l 



< Ctr 



3.6 X 10^'' T, 



25 rpl/2 , 



(4) 



finally, it must require a total amount of energy smaller 
than the total energy of a GRB: 



= 8.7 X 10-* Tg n tjr — < 10 
in 



52 



£^52 erg 



Conditions || and ^ yield the density constraint: 



n < 3.7 X lOTgS/flL 



46 



47r 



while conditions || and || give: 
n > 1.2 X 10^ T^^^ Lie E-^^ 



(5) 



(6) 



(7) 



The above conditions cannot clearly be satisfied simul- 
taneously if all the parameters are taken equal to their 
fiducial values. Since rjn and Sl/47r are both numbers 
smaller than unity, a change in the geometry does not help, 
making the constraint of Eq. [g] more stringent. Also the 
temperature cannot be changed significantly, especially as 
long as X-ray lines must be taken into account. 
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Fig. 3. Cartoon of the geometrical set-up for the GRB and 
thermal reprocessing material. The GRB is surrounded by 
a thin shell (or clumped thicker shell) of material, which 
is heated by GRB photons within the fireball solid angle. 

In conclusion, a thermal component that contributes 
to the early afterglow of a typical GRB cannot be emitted 
in the non-cooling regime. In fact, in order to reach a lumi- 
nosity large enough without involving a too large thermal 
energy, the plasma density must be so large to make the 
cooling time very short. As an example, R02 derived a 
cooling time of ~ 2 s for the afterglow of GRB 011211. 
We cannot therefore adopt the simple Eq. || but we must 
use the more complete Eq. |^. 

We therefore consider in the following two limiting 
cases: R/c{l — cos6) » icooi 3> thcat (flash heating) and 
-R/c(l — coa0) ^ thoat ^ ^cooi (steady heating). In prin- 
ciple also the case thoat ^ i?/c(l — cos 6') ^ tcooi may 
be interesting. It is however difficult to identify a heating 
source that can be active for a time comparable to the line 
emission time scale. Even if the emitting plasma would be 
heated by afterglow photons, the heating time could be at 
most a few per cent of the line emission time scale (Lazzati 
et al. 2002b). 

3.1. Flash heating 

Consider a portion of a spherical shell of plasma of radius 
R, temperature T and uniform density n covering a solid 
angle H. = 1— cos 9 (see, again. Fig. ^). Its emission integral 
can be written as EI = R^ rjB^ where -qn = AR/R 
is the ratio of the thickness of the shell to its radiusQ. 
Considering the values 7711,10 defined above (and shown in 



^ Note that if a clumpy medium is assumed instead of a 
spherical shell, riR is equal to the volume filling factor of the 
configuration that maximizes the emitted luminosity. 
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Fig. 4. The luminosity coefRcient Ch (see text) as a function of the plasma temperature and of metallicity 
considered elements. As in Fig. the dotted, dashed, solid, dot-dash and dot-dot-dot-dash lines are relative 
0.3, 1, 3 and 10 time solar metallicity, respectively. 



for the 
to 0.1, 



Fig.Jl|), the cloud will produce a line with luminosity^ (see 
Eq.B: 

Lii„c = 1.9 X lO^Vnci^Tgrr (8) 

where we have used the Thomson opacity tt = n Rrj^ ax- 
Eq. H coupled with the efficiencies shown in Fig. |l| implies 
that when we observe a line with luminosity iunc 5 x 
10'''* erg s^* (Ghisellini et al. 2002) a thermal model in 
the flash heating regime requires a radius larger than 
i? <L 2 X 10^^ cm. Analogous results are obtained by impos- 
ing a ten times larger continuum luminosity. Such a large 
radius implies that the line luminosity should remain con- 
stant for a time t = R(l — cos9)/c, to be compared with 
the observed variability time scales tvar ~ 2 x 10"* s ob- 
served in e.g. GRB 011211 (R02). This constraint can be 
satisfied if the heated plasma covers a very small solid an- 
gle, of opening angle 9 ~ 4.5°, at the limit but consistent 
with the smallest opening angles inferred for GRB jets 
from afterglow modelling (Frail et al. 2001; Panaitescu & 
Kumar 2002). In order to satisfy the cooling condition, 
however, the plasma must be dense (n ^ lO"'^ cm~^) and 
therefore confined in a very thin shell (or small clumps) 
with TjH ^ 5 X 10^^. In addition, if we assume that the 

^ The continuum luminosity can be easily derived by adopt- 
ing r;iine = 1. 



plasma is due to a spherical shell surrounding the GRB 
explosion site which is heated only in the small polar cap 
by afterglow photons, the shell mass would turn out to 
be M = 47r i?^ nip tt/<Jt ~ 600 Mq for the case explored 
above. Both the geometry factor rjn and the total mass re- 
quired seem to be rather extreme, even though the whole 
picture is, in this case, completely self consistent. 

3.2. Steady heating 

Let us now consider a plasma in which a source of heating 
is active for a time scale ihoat > 1-2 x 10*^ Tg/n = tcooi- 
In case of radiative heating, this time scale may be for 
example the GRB duration or even longer, if the afterglow 
photons can contribute to the heating. In steady heating, 
the line luminosity (continuum luminosity for T^iino = 1) 
can be written as: 

T o „ 1 o36 2 0^1/2 ^hcat ^hoat 

ilinc = 7.2 X lO'^^T^.Tg' ?7iino = Ch (9) 

where we have defined the parameter Ch in order to em- 
phasize the ratio of the heating time scale over the ge- 
ometry parameter. In Fig. ^ we show the highest possible 
values of the coefficient Ch as a function of temperature for 
five different values of metallicity (as in Fig. 0). The hues 
are plotted only when the equivalent width of the line is 



D. Lazzati: Thermal X-ray afterglows 



predicted to be larger than 100 eV, corresponding to the 
unshaded areas in Fig. |l|. 

To compare these results with afterglow data, consider 
a Fe Hne with L = lO^^ erg s^^ (GRB 991216; Piro et al. 
2000) or the S line in GRB 011211 (R02) with luminosity 
L — Ax 10^^ erg s^^. In both cases a ratio thcat/?/fl ^ 10^ 
is required, by comparison with the appropriate panel in 
Fig. ^. A typical heating time can be considered to be the 
fireball transit time, which is equivalent to the total du- 
ration of the prompt GRB emission. In both cases such 
duration is of the order of 100 s, leaving us with an ex- 
treme requirement of anisotropy rjj^ ^ 10^*". Alternatively, 
the heating may be provided radiatively by the absorp- 
tion of GRB and early afterglow photons. In that case, 
the duration of the heat supply may be of several per cent 
of the line emission time scale (Lazzati et al. 2002b). In 
that case the constraint would be relaxed to the (yet still 
challenging) value rjji ^ 10^^. The most convenient solu- 
tion is represented by an (unknown) heating mechanism 
acting for a time comparable to the line emission time 
scale itself. In this case, the constrain on the geometric 
factor rjfi would be similar to the one derived above for 
the flash heating case. Such a heating mechanism is how- 
ever presently obscure and will have to face the problem 
of stability discussed below (§4). 

In a steady heating scenario, therefore, thermal lines 
and continua can dominate the early afterglow emission, 
only in case of extreme geometric conditions, in which the 
emission is produced either in a sheet like shell or an ex- 
tremely clumpy medium. Interestingly, analogous extreme 
conditions were independently inferred for the environ- 
ment of GRB 000210 (Piro et al. 2002) in order to account 
for the lack of ionization features in the soft X-ray after- 
glow spectrum. Piro et al. (2002, and references therein) 
argue also that such conditions may be realized in giant 
molecular clouds. Stability considerations (see § 4), how- 
ever, show that more extreme conditions are required in 
this case. Differently from what derived in the flash heat- 
ing condition, in this case the radius of the shell does not 
have to be very large, and therefore the total mass required 
is not huge. 

4. Stability 

In the above section we have derived some geometrical 
constraints to the emitting plasma in order to reproduce 
the observed features and continua. Since the results envis- 
age particular geometric conditions, we here analyze their 
stability. First, since the heating energy is supposed to 
come from a central point (the GRB progenitor), one has 
to make sure that the emitting material is not accelerated 
to high velocities. 

Consider a shell of mass M absorbing energy from a 
relativistic outflow. If the material is radiative, it acquires 
a bulk velocity v — E/{Mc). A bulk velocity v ~ 0.1c 
was measured in GRB 011211 (R02) and GRB 991216 
(Piro et al. 2000). Requiring that our optically thin shell 
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is accelerated to a comparable or smaller speed implies a 
radius larger than: 

R > ^ = 3 X 10^3 cm (10) 

47r TOp rjline TT V 7/lino Tt Vq 

This is not a compelling limit, given the radii discussed 
above. 

Provided that the emitting medium is not accelerated 
to relativistic speeds by the energy input, we also want 
that the thin emitting shell (or blobs) do not expand in 
a time scale smaller than the emission one (which can 
be either the heating or cooling time scale). The shell (or 
blobs) was in fact in equilibrium with the ambient medium 
when it was cold. Now that its temperature is increased 
it will tend to expand under the effect of the increased 
internal pressure. If it expands at the speed of sound Cg, 
its density will be sizably modified in a time scale texp = 
Rrju/cs- It is therefore required that the expansion time 
is longer than the emission time scale. 

In the case of steady heating we obtain: 

> ^ = 10» t'^^ ^ > 10" t'/^ (11) 

Also in the steady heating case, therefore, the radius has 
to be large in order to allow for the production of ther- 
mal components in the early GRB X-ray afterglows. For 
such large radii, however, the density required to fulfill the 
steady heating conditions is n ^ lO^^/^hoat ^ 10"^^, several 
orders of magnitudes larger than what inferred by Piro et 
al. (2002). It seems therefore that the steady heating case 
requires more extreme conditions than the flash heating 
one (the same distance from the explosion site and lim- 
its on the total mass involved, but larger densities and 
smaller filling factors). Applying the same stability con- 
dition to a flash heating case, we are left with the more 
relaxed constraint 

TT > 0.1 (12) 

which is therefore not difficult to fulfill. These stability 
considerations suggest therefore that the only viable way 
to produce a sizable thermal afterglow component is by 
heating a small portion of a massive tt ~ 1 (either very 
thin or clumpy) shell of material located at a relatively 
large distance from the burst explosion site (Fig. ||). It 
should however be emphasized that these stability con- 
siderations cannot be applied if the heating is provided 
hydrodynamically, in such a way that the source of heat- 
ing is providing also the confining pressure. In this case 
also a steady heating scenario, with a much smaller shell, 
may be viable. 

5. dumpiness of GRB ambient medium 

The above discussion shows that the emission lines de- 
tected in several GRB afterglows can have a thermal ori- 
gin only if the GRB surroundings are extremely clumpy, 
with density and geometric contrasts of order of a hundred 
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thousand or more. We now discuss the reliability of these 
conditions in different GRB progenitor scenario. 

In the hypernova scenario, the GRB sets on simultane- 
ously to a supernova explosion, and the ambient medium 
is the result of the interaction of the pre-SN star with its 
surroundings. In particular, the nearby ambient will be 
dominated by the late stages of the mass ejection history 
of the star. These stages are known to be unsteady and 
clumps or shell-like structures can be envisaged. If the 
mass ejection is caused by radiative effects, however, ex- 
treme structure cannot be produced, since the stellar lumi- 
nosity varies on time scales comparable to the Kelvin time- 
scale, which is of the order of hundreds of years, yielding a 
thick shell. In the wind environment of SN1998bw, for ex- 
ample, Li & Chevalier (1999) found that inhomogeneities 
up to a factor of a few were present. A more appealing 
scenario, in this perspective, is the supranova model by 
Vietri and Stella (1998), where the GRB is supposed to 
explode several weeks to years after a supernova. The su- 
pernova explosion, having a much smaller time scale, can 
generate a more extreme geometry. The dumpiness of SN 
ejecta has been investigated by numerical simulations. It is 
found that Rayleigh- Taylor instabilities can produce high 
density clumps with angular scales ^ 1° (Bottcher et al. 
2002). Such structures have rjn = a/ R ^ 10~^ where a is 
the clump radius, much larger than the value required to 
produce sizable line emission from coUisional excitation. 
It is however possible, in the supranova scenario, that the 
remnant is illuminated by a super-Eddington relativistic 
wind (Vietri & Stella 1998, Konigl & Granot 2002) in the 
time span between the SN and GRB explosions. The inter- 
action of this wind with the ejected SN shell may increase 
the inhomogeneities originally present in the shell (Guetta 
& Granot 2003; Lazzati & Rees in preparation). 

6. Discussion 

In at least three XMM-Newton spectra of X-ray GRB af- 
terglows, the presence of a thermal component has been 
claimed (R02; Watson et al. 2002). We have studied the 
geometrical conditions that the thermal material must sat- 
isfy in order to contribute significantly to the early X-ray 
afterglow. We first concluded that the time scale during 
which the emission is observed must be set by the light 
crossing time of the emitting region. This imply that stan- 
dard free-free equations that relate the emission integral 
to the luminosity cannot be used and a more general equa- 
tion (Eq. |l|) has to be adopted. The implications of this 
are of great importance. In fact, the use of the standard 
formalism (Eq. |^) to compute the EI led R02 to underes- 
timate the emission integral for GRB 011211 by a factor 
~ 10^, and conclude that thermal emission from a shell 
with R ^ 10^^ cm can be a self consistent solution. We 
showed that, with the correct treatment of light-crossing 
effects, this is not the case. We then studied the possibility 
of (i) flash heating and (ii) steady heating. In both cases, 
we conclude that the thermal material must be extremely 
clumped in order to contribute significantly to a typical 



X-ray afterglow. In the steady heating case, the thermal 
material is continuously heated for a time longer than the 
cooling one, reaching a stationary state (albeit for a time 
smaller than the light-crossing one). We showed that in 
this case the material has to be so extremely clumped that 
the clumps are dissolved by the increased temperature in a 
very short time scale, so that an additional source of con- 
finement must be envisaged. In the case of flash heating, 
a self consistent solution can be found, requiring a less 
extreme (even if yet compelling) clumping. In this case, 
the thermal material must be located far from the source 
and therefore, in order to preserve the fast time variabil- 
ity, only a small portion of it has to be heated. Assuming 
that the same density is spread all over the GRB explo- 
sion site, we derive that the total mass involved is large 
(^ 600 Af©), larger than even a type-II supernova rem- 
nant. A moderate asymmetry of the remnant could how- 
ever ameliorate this requirement. It is interesting to note 
that, even though such clumpy and dense environments 
seem extreme, they were already inferred independently to 
account for the lack of ionization features in GRB 000210 
(Piro et al. 2002) and to explain the haigh energy emission 
in GRB 940217 (Katz 1994). 

To date there is not yet a compelling evidence that 
there are indeed thermal components in the early X-ray 
afterglow of GRBs. It is however intriguing then in three 
out of four afterglows observed by XMM-Newton a ther- 
mal model yields to a better fit than an absorbed power- 
law (R02, Watson et al. 2002a). In the case of GRB 011211 
(R02) is particularly difficult to explain the absence of an 
iron or nickel line, while reflection models seem to give a 
better explanation to the observed line ratios (Lazzati et 
al. 2002b). 

Nonetheless, should a thermal component be con- 
firmed in the X-ray afterglow of at least a sub-fraction 
of GRBs, it would be difficult to avoid the conclusion 
that (i) the GRB explosion site is surrounded by a mas- 
sive shell of matter, similar to a 1-2 years old supernova 
remnant (as expected in some versions of the Supranova 
model, Vietri & Stella 1998); (ii) that this remnant is ex- 
tremely clumped, suggesting that some precursor activity 
has been taking place. Whether such a dumpiness can help 
making the prompt emission in an external shock scenario 
(Dermer et al. 2000) heavily depends on the presence of 
clumping at smaller scales, a test that cannot be done 
with thermal emission from the clumps, since would con- 
tribute in any case a negligible flux to the early afterglow 
emission. 

Acknowledgements. I am grateful to G. Ghisellini, P. Kumar, 
M. J. Rees and F. Tavecchio for useful comments and discus- 
sions. I acknowledge financial support from the PPARC. 

References 

Arnaud, K. A., 1996, in ASP Conf. series, Astronomical Data 
Analysis Software and Systems V, eds. G. H. Jacoby & J 
Barnes (S. Francisco: ASP), 17 

Berger E. et al., 2000, ApJ, 545, 56 



D. Lazzati: Thermal X-ray afterglows 



7 



Bottcher M., Fryer C. L., Dermer C. D., 2002, ApJ, 567, 441 
Covino S., Lazzati D., Ghisellini G. et al., 1999, A&A, 348, LI 
Dermer C. D., Bottcher M., Chiang J., 2000, ApJ, 537, 255 
Frail D. A., Kulkarni S. R., Sari R., et al., 2001, ApJ, 562, L55 
Ghisellini G., Lazzati D., Rossi E., Rees M. J., 2002, A&A, 
389, L33 



Guetta P., G ranot J., 2003, MNRAS in press (astro 



ph/0208156) 



Konigl A., Granot J., 2002, ApJ, 574, 134 

Kumar P., N arayan R., 2002, ApJ submitted (astro 



ph/0205488D 



Lazzati D., Campana S, Ghisellini G., 1999, MNRAS, 304, L31 
Lazzati D., Covino S., Ghisellini G., 2002a, MNRAS, 330, 583 
Lazzati D., Ramirez-Ruiz E., Rees M. J., 2002, ApJ, 572, L57 
Li Z., Chevalier R. A., 1999, ApJ, 526, 716 
Liedahl D. A., Osterheld A. L., Goldstein W. H., 1995, ApJ, 
438, L115 

Meszaros P., Rees M. J., 1997, ApJ, 476, 232 

Mewe R., Gronenschild E. H. B. M., van den Oord G. H. J., 

1985, A&AS, 62, 197 
Panaitescu A., Kumar P., 2002, ApJ, 571, 779 
Piro L., Garmire G., Garcia M., et al., 2000, Science, 290, 955 
Piro L., Frail D. A., Gorosabel J., et al., 2002, ApJ, 577, 680 
Reeves J. N., Watson D., Osborne J. P., et al, 2002, Nature, 

416, 512 (R02) 

Rybicki G. B., Lightman A. P., 1979, '^Radiative processes in 
astrophysics", New York, Wiley-Interscience, 1979 

van Paradijs J., Kouveliotou C, Wijers R. A. M. J., 2000, 
ARA&A, 38, 379 

Vietri, M. & Stella, L. 1998, 507, L45 

Vietri M., Perola C, Piro L., Stella L., 1999, MNRAS, 308, 
L29 

Watson D., Reeves J. N., Osborne J., O'Brien P. T., Pounds 
K. A., Tedds J. A., Santos-Lleo M., Ehle M., 2002, A&A, 
393, LI 



